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Abstract During its Cycle I, the MAGIC telescope tar- 
geted about 250 hours several galactic sources sought 
to be, or detected previously by other experiments in 
the same energy domain, 7-ray emitters. This paper re- 
views some results of such MAGIC observations covering, 
among others, supernova remnants, the Galactic Cen- 
ter and microquasars. We will concentrate on the recent 
discovery at very high energy 7-rays of the microquasar 
LS I +61 303. 
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1 Introduction: The MAGIC Telescope 

MAGIC is a telescope for very high energy (VHE, E > 
50 — 100 GeV) 7-ray observation exploiting the Imag- 
ing Air Cherenkov (IAC) technique. It is located on the 
Roque de los Muchachos Observatory (28°45'30"N, 17° 
52' 48" W, 2250 m above see level) in La Palma (Spain). 
This kind of instrument images the Cherenkov light pro- 
duced in the particle cascade initiated by a 7-ray in the 
atmosphere. MAGIC incorporates a number of techno- 
logical improvements in its design and is currently the 
largest single-dish telescope (diameter 17 m) in this en- 
ergy band, yielding the lowest threshold (~50 GeV). It is 
equipped with a 576-pixel photomultiplier camera with 
a 3.5° field of view. MAGIC's sensitivity above 100 GeV 
is ~ 2.5% of the Crab nebula flux (the calibration stan- 
dard candle for IAC telescopes) in 50 hours of obser- 
vations. The energy resolution above 200 GeV is better 
than 30%. The angular resolution is ~ 0.1°, while source 

Javier Rico 

Institut de Ffsica d'Altes Energies 

Edifici Cn Universitat Autonoma de Barcelona 

08193 Bellaterra (Barcelona) Spain 

Tel.: +34 93 581 3551 

Fax: +34 93 581 1938 

E-mail: jrico@ifae.es 



localization in the sky is provided with a precision of 
~ 2'. MAGIC is also unique among IAC telescopes by its 
capability to operate under moderate illumination (i.e. 
moonlight and twilight). This allows to increase the duty 
cycle by a factor 1.5 and a better sampling of variable 
sources is possible. 

The physics program developed with the MAGIC tele- 
scope includes both, topics of fundamental physics and 
astrophysics. In this paper we present the results re- 
garding the observations of galactic targets. The results 
from extragalactic observations are presented elsewhere 
in these proceedings pp. 



2 Highlights of cycle I 

MAGIC's first observation cycle spanned the period from 
January 2005 to April 2006. About 1 /4 of the observation 
time (not counting that devoted to Crab nebula techni- 
cal observations) was devoted to galactic objects. The 
observations covered both, candidates and well estab- 
lished VHE 7-ray emitters, and included the following 
types of objects: supernova remnants (SNRs), pulsars, 
pulsar wind nebulae (PWN), microquasars (/iQSRs), the 
Galactic Center (GC), one unidentified TeV source and 
one cataclysmic variable. In this section we highlight the 
results obtained so far from such observations, and con- 
centrate, in sectional on the most interesting case of the 
microquasar LS I +61 303. 



2.1 The Crab nebula and pulsar 

The Crab nebula is a steady emitter at GeV and TeV 
energies, what makes it into an excellent calibration can- 
dle. The Crab nebula has been observed extensively in 
the past over a wide range of wavelengths, covering the 
radio, optical and x-ray bands, as well as high-energy 
regions up to nearly 100 TeV. Nevertheless, some of the 
relevant physics phenomena are expected to happen in 
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Fig. 1 Energy spectrum above 100 GeV from the Crab neb- 
ula measured by MAGIC in two different observation seasons. 
For comparison, the extrapolations down to 100 of HEGRA 
and Whipple measurements are shown. 



Table 1 Upper limits (UL) to the pulsed integral flux from 
three pulsars observed by MAGIC. The considered energy 
threshold (E t h), observation time (OT), assumed duty cycle 
(DC) and confidence level (CL) are also shown. 
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the VHE domain, namely the spectrum showing an in- 
verse Compton (IC) peak close to 100 GeV, a cut-off of 
the pulsed emission somewhere between 10 and 100 GeV, 
and the verification of the flux stability down to the per- 
cent level. The existing VHE 7-ray experimental data is 
very well described by electron acceleration followed by 
the IC scattering of photons generated by synchrotron 
radiation (synchrotron self Compton process). Probing 
the presence/absence of a small contribution of VHE 7's 
produced in hadronic interactions is a challenge for ex- 
perimenters. 

Along the first cycle of MAGIC 's regular observa- 
tions, a significant amount of time has been devoted to 
observe the Crab nebula, both for technical and astro- 
physical studies. The performance of the telescope has 
been experimentally evaluated and found in good agree- 
ment with the expectations and Monte Carlo simula- 
tions 00] . This has allowed us to perform routine analy- 
ses above 100 GeV, where the performance of our instru- 
ment is fully understood. On the other hand, a sample 
of 12 hours of selected data has been used to measure 
with high precision the spectrum down to ~100 GeV, as 
shown in figure ^ 21 • We have also carried out a search 
for pulsed 7-ray emission from Crab pulsar and two mil- 
lisecond pulsars @J|o]j albeit without positive result. The 
derived upper limits of the pulsed flux for the three ob- 
served pulsars are shown in Tabled 
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Fig. 2 Energy spectrum above 400 GeV of HESS J1813-178 
(top) and above 150 GeV of HESS J1834-087 (bottom), mea- 
sured by MAGIC. For comparison, also the spectra measured 
by HESS and that of the Crab nebula are shown. 



2.2 Supernova Remnants 

Shocks produced at supernova explosions are assumed 
to be the source of the galactic component of the cos- 
mic ray flux [SJ. The proof that this is the case could be 
provided by observations in the VHE domain. The ratio- 
nale is that the hadronic component of the cosmic rays 
-enhanced close to their source, i.e. the SNR- should pro- 
duce VHE 7-rays by the interaction with nearby dense 
molecular clouds. Although recent data seem to indicate 
that this is the case, it is difficult to disentangle the VHE 
component initiated by hadrons from that produced by 
Bremsstrahlung and IC processes by accelerated elec- 
trons. Therefore more data in the TeV regime together 
with multi-wavelength studies are needed to finally solve 
the long-standing puzzle of the origin of galactic cosmic 
rays. 

Within its program of observation of galactic sources, 
MAGIC has observed a number of supernova remnants. 
In particular, we are observing several of the brightest 
EGRET sources associated to SNRs, and the analysis of 
the acquired data is in progress. On the other hand, we 
have confirmed the VHE 7-ray emission from the SNRs 
HESS J1813-178 [7j and HESS J1834-087 (W41) 0. Our 
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Fig. 3 Sky map of 7-ray candidate events (background sub- 
tracted) in the direction of HESS J1834-087 for an energy 
threshold of about 250 GeV. Overlayed are 12 CO emission 
contours (black) [5| and contours of 90 cm VLA radio data 
(green) [Tig. 
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Fig. 4 Integral VHE 7-ray flux above 1 TeV as a function of 
time from the GC as measured by MAGIC during September 
2004 and May- July 2005. 



of the detected flux. Furthermore, there is no evidence 
for variability of the flux on hour/day time scales nor on 
a year scale, as shown in Figure 01 This suggests that 
the acceleration takes place in a steady object such as a 
SNR or a PWN, and not in the central black hole. 



results have confirmed SNRs as a well established popu- 
lation of VHE 7-ray emitters. The energy spectra mea- 
sured by MAGIC for these two sources are shown in Fig- 
ure [21 They are, both, well described by an unbroken 
power law and an intensity of about 10% of the Crab 
nebula flux. Furthermore, MAGIC has proven its capa- 
bility to study moderately extended sources by observing 
HESS J1834-087. The morphology of this object mea- 
sured by MAGIC is shown in Figure [3] Interestingly, the 
maximum of the VHE emission has been correlated with 
a maximum in the density of a nearby molecular cloud 
(shown in the figure by the contour lines of the 12 CO 
emission intensities). Although the mechanism respon- 
sible for the VHE radiation remains yet to be clarified, 
this is a hint that it could be produced by high energy 
hadrons interacting with the molecular cloud. 



2.3 Galactic Center 

We have also measured the VHE 7-ray flux from the 
GC The possibility to indirectly detect dark mat- 
ter through its annihilation into VHE 7-rays has risen 
the interest to observe this region during the last years. 
Our observations have confirmed a point-like 7-ray ex- 
cess whose location is spatially consistent with Sgr A* 
as well as Sgr A East. The energy spectrum of the de- 
tected emission is well described by an unbroken power 
law of photon index a — —2.2 and intensity about 10% 
of that of the Crab nebula flux at 1 TeV. This result 
disfavours dark matter annihilation as the main origin 



2.4 Other observation targets 

MAGIC observation program of galactic sources includes 
also other kinds of objects such as PWN, cataclysmic 
variables, unidentified sources and /iQSRs (both high 
and low mass) . The analyses regarding these sources are 
ongoing and will be reported in the near future. The one 
exception is the high-mass x-ray binary LS I +61 303, 
discussed in the section [3] 



3 The microquasar LS I +61 303 

3.1 Microquasars 

Microquasars are a subclass of stellar, x-ray binary sys- 
tems that display prominent radio emission, usually at- 
tributed to the existence of jets of relativistic particles. 
They are named after the similarities with active galactic 
nuclei (AGNs), since /iQSRs show the same three ingre- 
dients that make up radio-loud AGNs: a compact ob- 
ject, an accretion disc, and relativistic jets Hence, 
/iQSRs are galactic, scaled-down versions of an AGN, 
where instead of a super-massive black hole we deal with 
a compact object of just a few solar masses that accretes 
material from a donor star. The similarities with AGNs 
explain the large interest risen by /tQSRs. Crucial for 
our understanding of accreting systems is also that: (i) 
//QSRs are nearby objects, and (ii) they show very short 
timescale variability. Those reasons make these objects 
ideal laboratories for the study of the physical processes 
that govern how /LtQSRs and AGNs work. In particular, 
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the short timescale variability displayed by /iQSRs allows 
to see changes in the ongoing physical processes within 
typical time scales ranging from minutes to months, in 
contrast with the usual scales of years to observe such 
variability in AGNs. In addition, /xQSRs could measur- 
ably contribute to the density of galactic cosmic rays |13| . 

3.2 LS I +61 303 

One of the most studied /^QSR candidates is LS I +61 303. 
This system is composed of a compact object of unknown 
nature (neutron star or black hole) in a highly eccentric 
(e = 0.7) orbit around a Be star. The orbital period - 
with associated radio and x-ray JS] outbursts- is 
26.5 days and periastron passage is at phase 0.23 |16|. 
The phase and intensity of the radio outburst show a 
modulation of 4.6 years High-resolution radio imag- 
ing techniques have shown extended, radio-emitting struc- 
tures with angular extension of ^0.01 to ~0.1 arc-sec, in- 
terpreted within the framework of /iQSR scenario, where 
the radio emission is originated in a two-sided, prob- 
ably precessmg, relativistic jet [file = 0.6) 11|. How- 
ever, no solid evidence of the presence of an accretion 
disk (i.e. a thermal x-ray component) has been observed. 
LS I +61 303 is also one of the two /^QSR candidates po- 
sitionally coincident with EGRET 7-ray sources ^5] , and 
the only one located in the Northern Hemisphere -hence 
a suitable target for MAGIC. There are also hints of 
variability of the 7-ray flux |2D). However, the large un- 
certainty of the position of the EGRET source has not 
allowed an unambiguous association with LS I +61 303. 

3.3 MAGIC observations 

LS I +61 303 was observed in the VHE regime with 
MAGIC during 54 hours (after standard quality selec- 
tion, discarding bad weather data) between October 2005 
and March 2006 [22 • The data analysis was carried out 
using the standard MAGIC reconstruction and analysis 
software [7U51lll| . 

Figure |5] shows the reconstructed 7-ray map during 
two different observation periods, around periastron pas- 
sage and at higher (0.4-0.7) orbital phases. No signifi- 
cant excess in the number of 7-ray events is detected 
around periastron passage, whereas it shows up clearly 
(9.4cr statistical significance) at later orbital phases. The 
distribution of 7-ray excess is consistent with a point- 
like source and is located at (J2000): a = 2 h 40 m 34 s , 
8 = 61°15'25", with statistical and systematic uncertain- 
ties of ±0.4' and ±2', respectively, in agreement with the 
position of LS I +61 303. In the natural case in which the 
VHE emission is produced by the same object detected 
at EGRET energies, this result identifies a 7-ray source 
that resisted classification during the last three decades. 

Our measurements show that the VHE 7-ray emis- 
sion from LS I +61 303 is variable. The 7-ray flux above 



o 

m 

a 

62 00 



6130 



61 DO 



60'30 



62' 00 



61 30 



61 '00 



60 30 



2 l, 50 m 2 h 45 m 2 h 40 m 2 h 35 m 

RA 

Fig. 5 Smoothed maps of 7-ray excess events above 400 
GeV around LS I +61 303. (A) Observations over 15.5 hours 
corresponding to data around periastron (i.e., between or- 
bital phases 0.2 and 0.3). (B) Observations over 10.7 hours 
at orbital phase between 0.4 and 0.7. The number of events 
is normalized in both cases to 10.7 hours of observation. The 
position of the optical source LSI +61 303 (yellow cross) and 
the 95% confidence level contours for the EGRET sources 
3EG J0229+6151 and 3EG J0241+6103 (green contours) are 
also shown. The bottom right circle shows the size of the 
point spread function of MAGIC (la radius). From Albert el 
al. EU. 

400 GeV coming from the direction of LS I +61 303 (see 
Figure EJ) has a maximum corresponding to about 16% 
of the Crab nebula flux, and is detected at around phase 
0.6. The combined statistical significance of the 3 highest 
flux measurements is 8.7a, for an integrated observation 
time of 4.2 hours. The probability for the distribution of 
measured fluxes to be a statistical fluctuation of a con- 
stant flux (obtained from a x 2 fit of a constant function 
to the entire data sample) is 3 x 10~ 5 . The fact that 
the detections occur at similar orbital phases hints at a 
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Fig. 6 VHE 7-ray flux of LS I +61 303 as a function of 
orbital phase for the six observed orbital cycles (six upper 
panels, one point per observation night) and averaged for the 
entire observation time (bottom panel). Vertical error bars 
include la statistical error and 10% systematic uncertainty 
on day-to-day relative fluxes. Only data points wit h m ore 
than 2 a significance are shown, and 2a upper limits [22] are 
derived for the rest. The modified Julian date (MJD) cor- 
responding to orbital phase is indicated for every orbital 
cycle. From Albert et al. |21| . 



periodic nature of the VHE 7-ray emission. Contempo- 
raneous radio observations of LS I +61 303 were carried 
out at 15 GHz with the Ryle Telescope covering several 
orbital periods of the source. The peak of the radio out- 
bursts was at phase 0.7, i.e. between 1 and 3 days after 
the increase observed at VHE 7-rays flux (see Figured). 

The VHE spectrum derived from data between ^200 
GeV and ~4 TeV at orbital phases between 0.4 and 0.7 
(see Figure[H| is fitted reasonably well (x 2 /ndf — 6.6/5) 
by a power law function: 

dN/(dA/dt/dE) = (2.7 ± 0.4 ± 0.8) x 10~ 12 x 

^(-2.6±0.2±0.2 



cirr 2 s _1 TeV~ 



(1) 

where N is the number of 7-rays reaching Earth per unit 
area A, time t and energy E (expressed in TeV) . Errors 
quoted are statistical and systematic, respectively. This 
spectrum is consistent with that measured by EGRET 
for a spectral break between 10 and 100 GeV. The flux 
from LS I +61 303 above 200 GeV corresponds to an 
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Fig. 7 LS I +61 303 radio flux density at 15 GHz measured 
with the Ryle Telescope (green squares, right axis) and re- 
sults from the last orbital cycle observed by MAGIC (black 
dots, left axis), from 14 February to 8 March 2006. The MJD 
cor resp onding to orbital phase is indicated. From Albert et 
al. Im 
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Fig. 8 Differential energy spectrum of LS I +61 303 for 
energies between 200 GeV and 4 TeV and averaged for or- 
bital phases between 0.4 and 0.7, measured by MAGIC. The 
error bars show the la statistical uncertainty. The dashed, 
red line corresponds to the Crab nebula differential spectrum 
also measured by MAGIC. The solid, black line is a fit of a 
power law (also expressed mathemati call y in the inset) to the 
measured points. From Albert et al. |21|. 



isotropic luminosity of ~ 7 x 10 33 erg s _1 (assuming a 
distance to the system of 2 kpc [221) , of the same order 
of that of the similar object LS 5039 [21], and a factor 
~ 2 lower than the previous experimental upper limit (< 
8.8 x 10" 12 cm" 2 s^ 1 above 500 GeV) [25]. LS I +61 303 
displays more luminosity at GeV than at x-ray energies, 
a behavior shared also by LS 5039. 



3.4 Emission scenarios 

LS I +61 303 belongs, together with LS 5039 EU and 
PSR B1259-63 to a new class of objects, the so- 
called 7-ray binary systems, whose electromagnetic emis- 
sion extends up to the TeV domain. LS I +61 303 and 
LS 5039 are usually thought to be /iQSRs, since evi- 
dences of relativistic jets have been found at radio fre- 
quencies. In such scenario, the high energy emission is 
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produced by the shocks triggered at the relativistic jets 
|12| . However, no clear signal of the presence of an ac- 
cretion disk (in particular an spectral feature between 
~ 10 and ~ 100 keV due to the cut-off of the thermal 
emission) has been observed so far. Because of that, it 
has been alternatively proposed that relativistic parti- 
cles could be injected into the surrounding medium by 
the wind from a young pulsar |27| , which seems to be the 
case of PSR B1259-63. In the case of LS I +61 303 also 
the resemblance of the time variability and the radio/x- 
ray spectra with those of young pulsars seem to support 
such hypothesis. However, no pulsed emission has been 
detected from LS I +61 303. Therefore, the existing data 
in the radio/x/7-ray domain, cannot conclusively con- 
firm, or rule out, none of the two proposed scenarios. 
Thus, the question of whether the three known 7-ray bi- 
naries produce TeV emission by the same mechanism, 
and by which one, is currently object of an intense de- 
bate [2H|- Observation at VHE with MAGIC simulta- 
neously together with other instruments at other wave- 
length domains -in particular radio- will help elucidate 
the mechanism of the TeV emission and hence the nature 
of 7-ray binaries. 

The observation of jets has triggered the study of dif- 
ferent microquasar-based 7-ray emission models, some 
regarding hadronic mechanisms: relativistic protons in 
the jet interact with non-relativistic stellar wind ions, 
producing 7-rays via neutral pion decay [2H|; some re- 
garding leptonic mechanisms: IC scattering of relativistic 
electrons in the jet on stellar and/or synchrotron pho- 
tons 123]. 

The TeV flux measured by MAGIC has a maximum 
at phases 0.5-0.6 (see Figure EJ, overlapping with the x- 
ray outburst and the onset of the radio outburst. The 
timescale of the TeV variability constrains the emitting 
region to be smaller than a few 10 15 cm (or 0.1 arc-sec at 
2 kpc) , which is larger than the binary system and of the 
same order of the detected radio jet-like structures. The 
maximum flux is not detected at periastron, when the 
accretion rate is expected to be the largest. This result 
seems to favor the leptonic over the hadronic models, 
since the IC efficiency is likely to be higher than that of 
proton-proton collisions at the relatively large distances 
from the companion star at this orbital phase. Further, 
VHE 7-ray emission that peaks after periastron passage 
has been predicted considering electromagnetic cascad- 
ing within the binary system |31j . It is also possible to 
explain the orbital modulation of the VHE emission in 
LS I +61 303 by the geometrical effect of changes in the 
relative orientation of the stellar companion with respect 
to the compact object and jet as it impacts the position 
and depth of the 77 absorption trough [3*2j . The exist- 
ing data are, however, not conclusive to confirm or rule 
out any of the theoretical models existing in the litera- 
ture. Concerning energetics, a relativistic power of sev- 
eral 10 35 erg s _1 , extracted from accretion, could explain 



the non thermal luminosity of the source from radio to 
VHE 7-rays. 

LS I +61 303 is an excellent laboratory to study 
the VHE 7-ray emission and absorption processes taking 
place in massive x-ray binaries: the high eccentricity of 
the binary system provides very different physical condi- 
tions to be tested on timescales of less than one month. 
Future MAGIC observations will test both, the period- 
icity of the signal and its intra-night variability. 
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